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Abstract 

We study the eikonal model for the nuclear-induced breakup of Borromean 
nuclei, using ^^Li and ^He as examples. The full eikonal model is difficult 
to realize because of six-dimensional integrals, but a number of simplifying 
approximations are found to be accurate. The integrated diffractive and one- 
nucleon stripping cross sections are rather insensitive to the neutron-neutron 
correlation, but the two-nucleon stripping does show some dependence on the 
correlation. The distribution of excitation energy in the neutron-core final 
state in one-neutron stripping reactions is quite sensitive to the shell structure 
of the halo wave function. Experimental data favor models with comparable 
amounts of s- and p-wave in the ^^Li halo. 

I. INTRODUCTION 

Halo nuclei having a very weakly bound neutron pair (often referred to as Borromean 
nuclei) are interesting objects but they are difficult to study experimentally. Secondary 
interactions in radioactive beams have been an important tool, with Coulomb excitation 
providing quantitative data about the excitation properties [HH. Nuclear excitation is also 
important from an experimental point of view, but the theoretical interpretation of nuclear 
reaction cross sections deserves closer attention. In this work we attempt to make a link, 
as quantitative as possible, between the nuclear excitation observables and the fundamen- 
tal properties of a Borromean nucleus. The fact that correlations can play an important 
role makes this goal more difficult than for a nucleus with a single-nucleon halo. On the 
experimental side, we have been inspired by the work on ^^Li carried out at Ganil, NSCL, 
RIKEN and most recently at GSI. The extremely large Coulomb breakup cross section shows 
the halo character of the nucleus, but the details of its wave function have been controver- 
sial. Starting from the shell model, two of us @] constructed a wave function that fit many 
Coulomb excitation measurements [Q. It had a dominant shell configuration, as one 
expects from Hartree-Fock theory. However, several measurements (see for example ref. [^) 
and also the spectroscopy of the nearby nucleus ^^Be suggest a leading s^^g configuration in 
iiLi. 
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In principle, a nuclear-induced breakup gives independent information, and so it is de- 
sirable to calculate the various cross sections and compare with experiment. A recent ex- 
periment was carried out on a ^^C target at 280 MeV/n. At that energy it is justified to 
treat the target-projectile interaction in the sudden approximation, using the NN forward 
scattering amplitude for the interaction. Thus we may neglect the evolution of the wave 
function during the interaction time, provided we take the interaction from nucleon-nucleon 
scattering. The energy domain around 250 MeV has an additional advantage from a the- 
oretical point of view: the real part of the NN forward scattering amplitude goes through 
zero in this vicinity, so only the absorptive part of the interaction needs to be treated in the 
theory. 

The nuclear excitation of Borromean nuclei have been considered by a number of authors 



|0llT2|,|Tl],IT3[ . In treating the differential cross sections, it common to make a number of 



simplifying assumptions. We list them here: 

• Ground state wave function. Neutron-neutron correlations were neglected in ref. [§. 
We shall apply wave functions that have the full three-particle correlations. It turns 
out that differential cross sections are quite insensitive to these correlations, except 
the two-neutron stripping which does show an effect. Independent particle models can 
only describe pure configurations, so a mixture of s- and p-waves requires a correlated 
model. 

• Reaction model. In this work we use an eikonal model description of the nuclear 
reaction, improving on the black disk model of ref. 

• Neutron-core potential. It is important to include the final state neutron-core potential 
in calculating the energy or momentum spectra, as demonstrated in ref. [0]. Ref. [7] 
made use of simplified two-body wave functions which were based on a zero range 
neutron-core potential. We shall use a more realistic, finite range potential, both in 
the initial and final states. 

We shall investigate the validity of these as well as other approximations that are often 
made. Our main interest is the sensitivity of experiments to the properties of the halo 
nucleus. In a previous work |14] we developed models of the ^^Li ground state wave function 
with differing amounts of s-wave. One of our objectives is to see how well the amount of 
s-wave can be determined by the observables in a breakup reaction. The observables we 
consider are integrated cross sections for diffraction and one- and two-nucleon removal, and 
the differential cross section for the excitation energy in the ^Li + n final state when one 
neutron has been removed. 



II. REACTION MODEL 

The sudden approximation leads to the eikonal model for nucleus-nucleus interactions. 
In previous studies, we have applied the model to the nuclear-induced breakup of single- 



nucleon halo nuclei |T5|. Here we apply it to the breakup of a two-neutron halo nucleus. 
The effect of the interaction with the target to multiply the halo wave function by the 
profile functions S{Ri) for each particle, where i?j denotes the impact parameter of particle 
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i with respect to a target nucleus. The halo nucleus ^^Li has two neutrons and a ^Li core, 
requiring two types profile functions, Sn and Sc, associated with a neutron and the core, 
respectively. There are three integrated cross sections that leave the core intact, namely the 
diffractive, the one-nucleon stripping, and the two-nucleon stripping cross sections. These 
can be written 



(Tdif — J d^R, 



((1 - S,iR,)Sr,{Rl)Sr,{R2)f) - ((1 - S,(i?,)S„(i?i)S'„(i?2)))' 

{S'^iR,)SliR,)SliR,)) - {S,iR,)SniRi)S^iR2))' 



(la) 
(lb) 



ain-st = 2 J d^R,^{Sl{R,)Sl{R2){l - Sl{R^))), (2) 

a2n-st = J d'Rcm{S',{Rc){l - ^^(i?i))(l - 5'^(i?2))), (3) 

where Rem is the impact parameter of the halo nucleus with respect to the target nucleus, 
and (...) denotes a ground state expectation value. Our ground state wave function \l/o is 
expressed in terms of the relative neutron-core distances, ri and r2. An example of the 
needed expectation values is the one-nucleon stripping integral 



(5252(1 - SD) = I d'nd^r2\^o{ri,r2)\'S!{R,)Sl{R, + r2x)(l - Sl{R, + 



(4) 



The integrations are here performed for fixed Rem so one should use Rc = Rem — (^i± + 
1^21.) /{Ac + 2), where Ac is the mass number of the core nucleus. 

The six- dimensional integration in eq. (^) is very time consuming to carry out unless 
some simplification are made in the wave function or in the profile functions. We shall 
consider two simplifying approximations. The first is the no-recoil limit in which the impact 
parameter of the core, Rc, is assumed to coincide with the impact parameter Rem of the 
halo nucleus. The core profile function Sc can then be taken outside the expectation value. 
In addition, the integrations over ri and r2 become independent and the six- dimensional 
integral separates into a product of two three-dimensional integrals.^ Another simplifying 
approximation is the transparent limit, defined here by setting the factor S^{R2) equal to 
one inside the expectation value of eq. (0), thus neglecting the absorption of the second 
neutron. These two assumptions yield the cross section 

a,n-st,trans = 2 j d^ R Sl{R){l - Sl{R + T^^)) . (5) 

Note that this cross section is identical to the sum of the one- and two times the two-neutron 
stripping cross section, 

^ln~st,trans = + 2(J2„_st. (6) 



As discussed later on, the no-recoil limit differs from the exact calculation only in the case of 
diffraction. 
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We will see later that the two-neutron stripping cross section is rather small, so the trans- 
parent limit is a good approximation for this cross section. 

Of course, much more information about the halo is contained in differential final state 
distributions. The diffractive cross section has three particles in the final state, but that 
distribution is beyond what we can calculate, requiring three-particle continuum wave func- 
tions for many partial waves. The one-neutron stripping leaves two particles in the final 
state, and the differential cross section for that state is amenable to computation. The 
expression for the momentum distribution associated with the relative motion of the two 
surviving particles is 



da 



St 



d^k 



2 / d^R, 1 - SliRi) / dW M{R,,r2,,k) 



(7) 



where r2c is the center-of-mass coordinate of the remaining neutron-core system with respect 
to the stripped neutron; the associated impact parameter with respect to the target nucleus 
is denoted by i?2c, R2c = Ri + r2cj_- The amplitude M is given by 



M = J d\2 rc{k,r2)S,{R,)Sn{R2)^o{rur2). 



(8) 



Here ipc{k,r2) is the continuum wave function of the surviving neutron-core system, nor- 
malized to a plane wave at infinity. The coordinates Rc, R2 and ri are expressed in terms 
of the integration variables as Rc = R2C — T2±/(^c + 1), R2 = Ric + T2±Ac/{Ac + 1), and 
n = -r2c + r2/{Ac + 1). 

The numerical calculation of eq. (H) is rather difficult because of the form of the ground 
state wave function that we apply (see the next section). A major simplification is achieved 
by adopting the approximation ri = —r2c in the ground state wave function. The amplitude 
is then given by 



M' = J d\2 rAk,r2)Sc{Rc)Sn{R2)-^o{-r2c,r2). 



(9) 



An even simpler approximation is to also ignore the recoil correction in the argument of the 
core profile function, i. e. setting R^ = R2C, and to use the transparent limit for the second 
neutron, i. e. setting S'„(i?2) = 1- We shall refer to these approximations as the no-recoil, 
transparent limit, where the amplitude reduces to 



M" = Sc{R2c) / d\2 rc{k,r2)^o{-r2c,r2). 



(10) 



This approximation is used in ref. M and |T3. We will discuss the validity of the various 



approximations in Sec. VII below. 



III. THE THREE-BODY WAVE FUNCTION 



In ref. we constructed several three-body models of ^^Li. The models are based on 
Hamiltonians that all reproduce the empirical neutron-neutron scattering length and all 
have a binding energy of the three-body system close to the empirical value of 295 ± 35 
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keV 10 . The single-particle potentials and the density-dependence of the neutron-neutron 
interaction are varied to produce different probabilities of s- and p-waves in the different 
models. We will specifically examine the observables for models having 4.5%, 23%, and 50% 
s-wave. The first model is similar to the one used in ref. 0. The second and third models 
are constructed with a deeper neutron-core potential for even parity, single-particle states 
to increase the s-wave component in the ground wave function. Details of the procedure and 
the first two models are given in ref. WM. The wave functions are calculated in the form of 



single-particle states u{r) and amplitudes a as 

^oiri,f2) =J2J2'^ljnn'Uljnin)uijn'{r2) (11) 
Ij n,n' 

where r*i and r2 are neutron-core separation vectors. The indices [IsY label the single- 
particle, spin-angle wave functions which are coupled to zero total angular momentum as 
indicated by the superscript on the bracket. The indices n,n' label the radial quantum 
numbers of the single-particle basis states. These states are discretized by putting the 
system into a spherical box of large radius (typically 40 fm). 

In Table I we give some characteristics of the Hamiltonians and the resulting ground state 
wave-functions. A useful quantity is the s-wave scattering length. Within the constraints of 
our three-body model, we can only produce wave functions that are predominantly s-wave 
by using neutron-core potentials that produce extremely large s-wave scattering lengths. 
The parameters of the potential for the odd-parity states are fixed by the position of the 
Pi/2 resonance, which we assume to be at Er = 540 keV as suggested by measurements 



1% . Besides the s-wave probability, these wave functions show significant differences. The 
single-particle densities of the three models are shown in Fig. |I|. It may be seen that the 
halo is more extended the larger the s-wave probability. This is also apparent from the mean 
square neutron radii computed in Table |I|. 

The integrated dipole strength for Coulomb breakup is proportional to the mean square 
radius of the two-neutron center of mass, given in the last column of Table I. The value 
obtained with the p89 wave function is consistent with the experimental Coulomb breakup, 
but the s50 value is much too high. Thus, we cannot regard that wave function as realistic. 

Another important property of the wave function is the correlation between the two 
neutrons. The integrated dipole strength is proportional to the dineutron-core mean square 
radius, which in turn depends on the matrix element of fi- r2, as shown in [Q. In that work 
it was found that the correlation increased the dipole strength by 43 %. This enhancement 
does not depend very much on the model; with the present wave functions it is in the range 
of 30-40%, which may be extracted from the mean square radii computed in Table |. 

We also constructed s- and p-wave independent-particle models for comparison purposes. 
In these models, the single-particle potential is adjusted to match the exponential fall-off 
of the single-particle density that is obtained with the three-body model s23 mentioned in 
Table |. We see from Fig. |I] that the pure s-wave model gives a much more diffuse halo than 
the p-wave model or the correlated models. 
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IV. PROFILE FUNCTIONS 



We now specify the profile functions Sn and Sc that we use for our cross section calcu- 
lations. The neutron profile function Sn in the eikonal approximation is expressed in terms 
of the density of the target pt and the nucleon-nucleon cross section as 



S{b) = exp 



0"r, 



(12) 



We model the density of the C target with the harmonic oscillator fit to the charge density 



of ref. 118 



p(r) = po 1 + a(r/a)^ 



,-(r/a)2 



(13) 



with a = 1.687 fm and a = 1.067. The nucleon-nucleon cross section is taken from ref. [lU 
it is 29.2 mb at 280 MeV beam energy. 

The reliability of the model can be checked against the nucleon-carbon cross sections. 
The predicted reaction and elastic cross sections in the eikonal model are 



1 - slib) 

1 - Snib) 



(14) 
(15) 



These are compared with experiment in Fig. |^. The nucleon-carbon reaction cross section is 
taken from the proton cross section data tabulated in ref. ||2^, quoting ref. 0]. The total 
cross section for nucleon-carbon scattering is taken from the neutron measurements of ref. 
2^ . The experimental elastic cross section is deduced from the difference between total and 



reaction cross section. The agreement between our parameterization of Sn and experiment 
is close enough that we will not attempt to adjust the profile function to make a better 
fit. In the next section we will discuss how the cross sections in halo nuclei depend on the 
nucleon-target cross sections. 

The core-target profile function requires the convolution of both densities, 



Sc{b) = exp 



J dxdy J dzpc[\l {x — 6)^ + + -2^) j dz' pt(^\J x"^ + y"^ + z'"^ 



For the density of ^Li, we note that it has the same number of neutrons as ^^C and we will 
accordingly take the same parameters for the neutrons. The proton density does not have as 
many particles in the p shell, and we apply the pure harmonic oscillator model to determine 
a (= 1/3), and keep a the same as in ^^C. The resulting ^Li density is parameterized as 
in eq. (|T3|) with a = 1.687 and a = 0.726. This model gives an rms charge radius of 2.28 
fm, slightly smaller than the empirical charge radius of ''Li which is 2.39 fm. However, the 
predicted cross section at 800 MeV/n is 840 mb (assuming (T„„ = 40 mb), just 5% larger 

The cross section at 280 MeV/n 



23 



than the measured cross section of 796 ±6 mb from ref. 
has been measured for the mirror nucleus by Blank et al. 



They find a cross section 



of 812 ± 34 mb to be compared with 796 mb obtained by our model. 
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For the ^He core of ^He we use a 3-parameter- Fermi density function , 



Pc{r) = (1 + wr"^ /c^)/ [1 + exp((r - c)/z)\ , 

with w = 0.517, c = 0.964 fm, and z = 0.322 fm. At 800 MeV/n we find a total cross section 
of 546 mb (486 at 280 MeV/n) again comparable with the experimental result of 503 ±5 mb 



V. INTEGRATED CROSS SECTIONS 

In this section we examine the integrated cross sections and compare with experimental 
data. It is important to understand the dependence of the cross sections on the interaction 
and on the properties of the halo, and to this end we first discuss some estimates and bounds 
on the cross sections. 



A. Cross section bounds 

Let us first consider the (unrealistic) limit where the core-target interaction is ignored^. 
For a single-nucleon halo the stripping cross section would then be identical to the reaction 
cross section are between the nucleon and the target. In the case of a two-nucleon halo, 
the stripping cross section is just doubled. Split into In- and 2n-stripping components, the 
relation is 



cr 





In— st 



+ '2<-st = 2fT.e. (16) 



The symbol o"° is a reminder that the core shadowing is neglected, i.e. the factor Sc is set 
equal to one. Eq. (^) is illustrated in Fig. 0, showing the comparison of the left- and 
right-hand sides of the equation for the case of ^^Li breakup on a ^^C target. We used 
the s23 model to evaluate the unshadowed cross section. The relative amounts of In- and 
2n-stripping depends of course on the wave function and details of the interaction; in the 
case considered here, the In-stripping cross section is an order of magnitude larger than the 
2n stripping. 

The diffraction cross section is much more difficult to bound or estimate without full 
calculation of the integrals. In the case of a one-nucleon halo, a bound can be obtained 



by dropping the second term in the equation analogous to eq. ([Tal). The first term is just 



the elastic nucleon-target cross section, so the bound is a^^j < a^i- For Be at 800 MeV/n 



and using the model of |T5| we get a^^r = 68 mb, smaller than the bound aei = 85 mb. 



This shows that the neglected term is not necessarily small. For the two-nucleon diffraction 



formula eq. (|1^) the first term does not reduce to the elastic scattering, so no strict bound 



can be obtained this way. 



^This is commonly referred to as the transparent limit but we have reserved that concept for the 
transparency of the second neutron in a In-stripping reaction, c. f. Sec. 
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Instead, we shall analyze diffractive cross section qualitatively using a grey-disk model 
for the nucleon-target interaction. Thus, we take the nucleon profile function to be of the 
form Sn{b) = 1 — t9{bQ — b), < t < 1. We also need to assume that the two nucleons are 
uncorrelated in the halo wave function. The diffraction cross section may then be expressed 

as 

alj = 2a Jl - (1 + 2t - tV2)a(pt) + 2ta^(p?) - ^^a'{p',)] , (17) 



where (p^) is the ri-th moment of the transverse nucleon density in the halo^ and a = nb^ 



2 

2 


The leading term is twice the nucleon-target elastic cross section and the corrections are 
controlled by the parameter a{pt). The coefficient of the first correction is negative, so the 
first term gives a bound that is valid for large halos and small targets, 

^If < 2^d. (18) 

In the case of a ^^C target, the experimental elastic and reaction cross sections may be fit 
with a = 30 fm^ and t = 0.5. The transverse halo density for the ^^Li wave functions has 
the order of magnitude (pt) ~ 1/100 fm~^ and the second term makes about a factor of two 
correction; the higher terms are less important. The actual numbers for our model of the 

.0 

dif 

factor 2/3, as expected from the above analysis. 



^^Li-^^C reaction are shown in Fig. R|. We find a^^f = 75 mb, reduced from 2aei by about a 



B. Core shadowing 

In this section we examine the effect of the core shadowing on the cross sections, and 
use again the model s23 for the ground state wave function. Fig. |^ shows on the right the 
shadowing effect of the carbon target in the ^^Li breakup reaction. The In-stripping cross 
section is reduced to 43% of the unshadowed value cr'^n-st- The shadowing factor for the 
diffractive cross section is very similar (44%) to that for the In-stripping. 

The shadowing factor for two-nucleon stripping is much stronger than for the other 
processes; it reduces the cross section to 20% of the unshadowed value. The difference may be 
understood qualitatively as follows. The one-nucleon stripping and the diffractive excitation 
require avoiding an absorptive interaction with a least one of the halo nucleons, favoring 
moderately large impact parameters. On the other hand, the two-nucleon absorption has no 
such restriction and would be concentrated entirely at small impact parameters but for the 
presence of Sc- The different dependences on impact parameter is shown in Fig. ^. Here we 
see that the 2n-stripping probability is more concentrated at small impact parameter than 
the In-stripping and the diffractive probabilities, which are very similar to each other. 

The shadowing factor varies, of course, with target size. This dependence is illustrated 
in Fig. ^, where the target densities were taken from We see that the shadowing 

changes by a factor of 2 both for diffraction and In-stripping, going from a ^He target to a 
heavy target, and by a factor of 4 for the 2n-stripping. 



More precisely, it is the n-th moment of an averaged transverse density, the averaging being over 
the shape of the nucleon profile function. 
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C. Wave function sensitivity 



We next consider the sensitivity of the cross sections to properties of the halo wave 
functions. The various cross sections for different models are given in Table ||. For the single- 
nucleon stripping cross section, the shadowing factor varies, depending on how extended the 
single-particle density is. From Table |, we see that the mean square radius of the halo 
increases as the s-wave probability increases. Thus we expect less shadowing and a larger 
cross section for the models with a larger s-wave. This is indeed born out by the numbers 
in Table |I|. For the two-nucleon stripping, the correlation between the neutrons should be 
important as well, as they must both interact with the target. Indeed, we see from Table |I| 
that the two-nucleon cross sections doubles going from an uncorrelated p-wave model to the 
model with correlations and p-wave dominance. This may be compared with the effect of the 
correlations on the dipole transition strength which, as was mentioned in the last section, 
gives only a 30-40% enhancement. 



D. He cross sections 



Here we report corresponding cross sections for the breakup of ^He on a ^^C target, 
using the ®He wave function from ref. line 5 of Table II. The cross sections for two 
different beam energies are shown in Table |I|. ^He is more tightly bound than ^^Li, so 



the halo density does not extend out as far. Another difference is that He has a dominant 
P3/2 shell configuration which allows a stronger spatial correlation; pure Si/2 or pure pi/2 
configurations, on the other hand, have uncorrelated densities. The larger correlation implies 
that the 2n-stripping will be relatively stronger. This is indeed seen to be the case in 

the 2n-stripping cross sections is about a factor of two larger for ^He than for 

Otherwise, the cross sections are about the same as for ^^Li. The shadowing factors 
are similar, due to balancing features of a smaller core and a less extended halo. 



Table |H 
"Li 



VI. SHAPE OF THE STRIPPING SPECTRUM 

In this section we discuss analytic forms for the shape of the spectrum in the neutron-core 
system produced by the In-stripping reaction. The standard parameterization of a peaked 
distribution by the Breit-Wigner function is not justified at energies close to zero, or for the 
overlaps with extended wave functions. We shall propose parameterizations that take into 
account the threshold behavior and the halo character of the initial state. Throughout this 
section we make use of the no- recoil, transparent limit defined in eq. ([T0|). 

To treat the one-neutron removal from a Borromean nucleus, we simply take the overlap 
of the initial ground state wave-function with the continuum final state of the neutron-core 
system (cf. eq. (pUD), fixing the position ri = —r2c of the stripped neutron. The stripping 
model assumes that the process is incoherent in ri. Thus we consider matrix elements of 
the form 

M"{n,k) = J ^fc(r2)^o(ri,r2) (19) 
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and a probability distribution of the form 



\M"{n,k)\^dnk. 

Here dn^ ~ k'^dk for a differential momentum distribution and drik ~ kdE for a differential 
distribution in excitation energy of the neutron-core system. 

Different partial waves of the continuum wave function are incoherent if we integrate 
over the direction of the decay distribution. At low energies the s- and p- waves will be most 
important, and we now discuss their functional behavior. 

A. s-Wave distribution 

The s-wave distribution can be described analytically in the limit where the wave func- 
tions are dominated by their asymptotic behavioiQ. The continuum s-wave is then given 
by 

sin(A;r -|- 6) 



kr 



Here 6 is the s-wave phase shift; the scattering length is the linear coefficient in the expansion 

S = -ak + 0{k^). 

The two-particle initial state wave function has no exact analytic limit, but the general 
exponential fall-off at large distances suggests the approximation 



^(r2,r) ^ /(r2)- 



Then the integral in eq. ([T9|) can be carried out to give 

^„ cos(/ca) — a/ksm.{ka) 
k"^ + 

and therefore the cross section for the s-wave is given by 



da 
dE 



- 1 


2 r 


[a^ + k^\ 





a 



cos(ka) — — sin(A;a) 
k 



(20) 



This depends on the initial state potential through the fall-off parameter a and on the final 
state potential through the scattering length a. If the two potentials are the same, then the 
orthogonality of initial and final states requires the matrix element to vanish. This comes 
about in eq. (20) to leading order in k by the well-known relation between binding energy 



and scattering length [28|. 



^The formula derived gives an excellent description of the magnetic photo-disintegration of the 



deuteron at low energy |27|| . 
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The s-wave energy distribution for ^^Li is shown in Fig. ^ Here we have fitted both 
parameters a and a to give the best agreement with the calculated curve, which was obtained 
in the no-recoil, transparent limit (cf. eq. 10) using the model s23. In principal a is given 
by the binding energy as a ~ \/2mEB- Our fit has a = 24.5 MeV/c. The corresponding 
binding energy is 0.32 MeV, almost equal to the binding energy of the independent particle 
model. Also our fitted value a = —4.7 fm is very close to the actual scattering length of 
—5.6 fm. 

The distribution in Fig. |^ peaks at very low energies; the peak position is close to the 
momentum k = a/2 for a fairly wide range of scattering lengths a between —1/a < a < 1/a. 
This corresponds to an energy peak Epeak at 

Epeak = Eb/4:. 

With our theoretical fit, Ep^ak ~ 0.08 MeV. It should also be mentioned that for models 
with very large scattering lengths, such as the s50 model, the scattering length sets the 
momentum scale and the predicted peak is lower in energy. 

The corresponding s-wave distribution for ^He is shown in Fig. ^ Here the best scattering 
length parameters is a = 1.6 fm, to be compared with the actual scattering length of a = 2.4 
fm associated with the ^He potential. The best fit value of a is a = 55 MeV/c; this may be 
compared with the binding energy estimate \/2rnEB = 43 MeV/c. 



B. p-Wave distribution 

For the p-wave, measurements of the ^^B('^Li,^B)^°Li reaction have suggested the exis- 
tence of a resonance at about 540 keV Jl^. In our recent study of the ^^Li wave function we 



used this data to fix the p-wave potential for the neutron-core system. In this section we 
wish to establish a simple function to represent the distributions that we calculate. After 
trying different functional forms, we found that one could get acceptable fits with the Breit- 
Wigner resonance form but using an energy-dependent width. The threshold behavior of a 
p-wave resonance requires a width depending on energy as F ~ E^^"^. However, the width 
cannot continue to grow as the 3/2 power at energies above the resonance. We shall account 



for this by using the form of the p-wave width obtained in potential scattering p9 

The Breit-Wigner function for the decay of a resonance is then given by 

— = A ^ (22) 

dE {E -ErY + T^/A ^ ' 

There are four parameters here, namely the resonance energy Er, the width on resonance 
Fq, the overall strength A, and a cutoff parameter c. One might think that Fq and c 
could be determined by the radius of the potential forming the resonance, but the fact that 
we are considering a transition from a halo state implies that the length scales are larger 
than the nuclear radius. We shall treat them as adjustable parameters. Fig. ^ shows a fit 
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with parameter values Fq = 1 MeV, = 0.51 MeV, and c = 1.7. When we make an 
unconstrained fit, the parameters Fq and c acquire large values, showing that the function 
F is close to the E^^'^ dependence, except at extremely low energies. The deviation from 
the 3/2 power law is quite understandable; although the single-particle wave function is in a 
p-wave, the fact that the initial state is a halo means that the extension of the wave function 
rather than the radius of the barrier sets the scale for the threshold region. 

In Fig. we show similar a comparison for ^He stripping. In this case, the peak of the 
p-wave distribution is located at 0.83 MeV which corresponds quite well to the resonance 
energy of the P3/2 scattering state (0.89 MeV). Nevertheless, the best fit again favors large 
values of c, indicating that the 3/2 power law for the width is only valid very close to the 
threshold. 

VII. STRIPPING DISTRIBUTIONS 

In this section we treat the full stripping cross section as defined by eqs. (|^) and @. 
There are a number of questions one can ask: i) How reliable is the transparent limit to 
deduce the shape of the differential cross section? Do we get the right shape for the s and 
p-wave and also the right ratio between them? ii) Can a decomposition of the differential 
cross section into s- and p-wave components be used to infer the s/p ratio in the ground 
state wave function? 

The energy distribution for the full calculation is compared in Fig. ^ to the no-recoil, 
transparent limit for ^^Li stripping, using the s23 ground state wave function. We see that 
the effect of the neutron shadowing is to reduce the cross section without affecting the 
energy distribution. Thus we can use the transparent limit with confidence in describing 
these distributions. 

The In-stripping calculation involves integrations of the probability over the impact 
parameters of the absorbed neutron and the center of mass of the second neutron with 
respect to the core. These integrals are incoherent, and can be performed on a coarser 
mesh than the wave function integral over the internal coordinate. It is of interest to see 
how much the shape depends on the impact parameters, to make this integration as coarse 
as possible. Fig. ^ shows the distribution with the s23 wave function, comparing the full 
calculation with the distributions at typical impact parameters. Here we fix the value of the 
transverse radii i?2c and -Ri and integrate over the angle between them, and also over the 
longitudinal component of r2c- The curves have been normalized to the same area. We see 
that the ratio of cross sections at low energy to that at the p-wave peak varies by the order 
of 20%. Thus the impact parameter integration mesh cannot be coarser than 1 fm or so. 

Next we consider the relation of the s- and p-wave probabilities in the wave function 
to the corresponding probabilities in the stripping distributions. In principle, the stripping 
distributions for s- and p-waves are sufficiently dissimilar that one should be able to extract 
their relative probabilities from the data. In Fig. ^ we show the percentage of these waves 
as would be found in the energy distributions for the full calculation and also in the no-recoil, 
transparent limit, using the s23 wave function. As found before, the transparent limit is 
very accurate for purposes of extracting these ratios. As could be anticipated from Fig. |^, 
it is necessary to do a complete integration over all impact parameters to determine the 
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theoretical ratios. Finally, the amount of s-wave in the final state stripping distribution is 
systematically larger than in the initial ground state wave function. This is due to the larger 
extension of the s-wave, and the resulting decreased shadowing by the core. In Fig. ITT] we 



compare the s- and p-wave probabilities of the three correlated models with the final state 
probabilities for In-stripping. In all cases the s-wave probability is somewhat larger in the 
final state than in the wave function, for the reason given above. 



VIII. THE NO-RECOIL APPROXIMATION 



In this section we examine the no-recoil approximation, and will find that it is very 
accurate for the ^^Li reaction. The no-recoil approximation is in fact exact for the integrated 
stripping reactions, because the outer integration over Rem can be changed to an integration 
over the core coordinate by a simple change of variable. This separation cannot be made for 
the diffractive cross section, so there will be some effect of including core recoil. We discuss it 
immediately below. Following that, we discuss the energy distribution of the single-particle 
stripping, which is also subject to recoil corrections. 



A. Diffraction 



Diffraction arises from the fluctuations in the expectations of powers of the proflle oper- 
ators (see eq. (lb)), and will only be signiflcant at impact parameters where the operators 
vary. There are two effects that make the diffractive scattering dependent on the distinction 
between the center-of-mass coordinate and the core coordinate. The flrst is that there are 
no fluctuations associated with the core proflle function in the no-recoil limit, and includ- 
ing those fluctuations might be expected to increase the diffractive cross section. However, 
there is another effect that goes in the opposite direction. That is that the spatial extension 
of the neutrons is smaller measured with respect to the center of mass than with respect 
to the core. One might expect a smaller fluctuation of Sn if the neutron wave functions 
are more conflned. Another way to look at the question is with eq. (0). The first term 
in this formula can be evaluated in either coordinate system. The effect of the no-recoil 
approximation is thus confined to the evaluation of the second term. We have evaluated the 
integral numerically for the uncorrelated p-wave model, and we find that the effect of recoil 
is to reduce the diffractive cross section by about 20%. Since the corrections are expected 
to scale as 1/A, one should not ignore them for nuclei lighter than ^^Li. 



B. Stripping distributions 



As shown in Sec. |T|we include the main recoil effect (the difference in the profile functions 
between Rc and i?2c) in our calculation. The only approximation is then to replace the 
position ri of the stripped neutron by — r2c, cf. eq. (|]). As this neutron is absorbed, whereas 
the other neutron and the core must survive, the major contribution to the spectrum comes 
from a region where ri is not small. Therefore we do not expect to get a big effect from this. 
Nevertheless, we have also calculated the exact spectrum from eq. (^ in the independent 
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j9-wave model for fixed impact parameters i?2c and Ri. The overall relative deviation of the 
two spectra is smaller than 3% with deviation near the peak being only of the order of 1/2%. 
The recoil effect on the spectrum therefore seems to be negligible. 



IX. COMPARISON WITH EXPERIMENT 

A. Integrated cross sections 

In Table |I| we compare the integrated cross sections with the experimental data of ref. 
0, Table 3. The yields of events having zero, one, or two neutrons in coincidence with 
the ^Li fragment give respectively the cross sections labeled 2n-stripping, In-stripping, and 
diffraction. The first thing to note in the comparison with theory is that the total two- 
neutron removal cross section measured, 280 mb, is much larger than the eikonal model 
predicts. There are also data on the two-neutron removal cross section at 800 MeV/n 
pO| . Applying the eikonal model[| at this energy gives a cross section of 187 mb (s23 wave 



function), only 15% lower than the experimental value of 220 ± 10 mb. At the lower energy 
the s23 model gives 179 mb. This is close to the theoretical value at the higher energy, 
which is certainly to be expected in view of the mild change in the nucleon-nucleon cross 
section between the two energies. The direction of the change in both the nucleon-nucleon 
and the eikonal removal cross section is a decrease at the lower energy. In contrast the 
experimental value is larger at the lower energy. The theoretical two-neutron removal cross 
section behaves the same way in the case of ^He, as may be seen in the bottom line of Table 
III. Here also there is fair agreement with the experimental value at the higher energy. 

Let us now turn to the individual components. The diffractive and In-stripping cross 
sections are within experimental error of the most extreme wave function, s50, but the 2n- 
stripping cross section far exceeds any of the models. In principle, additional contributions to 
the cross section could come from processes outside the scope of the eikonal model. The flux 
that is absorbed in the eikonal can reappear, e.g. by multistep or rearrangement processes. 
This might first show up in cross sections that are very small in the eikonal model-such as 
the 2n-stripping. On the experimental side, the definition of the cross sections requires some 
momentum cut on the neutrons to be considered as part of the projectile fragmentation. In 
ref. [6], the detector acceptance allowed neutrons with transverse momenta of up to 60 
MeV/c to be included; the authors also quoted cross sections assuming that 20% of all 
fragmentation neutrons were outside that momentum range. The deduced cross sections are 
given in the last column of Table III, labeled "modified data". 

As may be seen from the table, adding the detector acceptance correction explains the 
small 2n-stripping cross section, but only at the expense of greater disagreement with theory 
for the other two cross sections. Thus, for the diffractive cross section, we obtain only 
between 33% and 51% of the experimental value. We believe that it not possible to explain 
this discrepancy within the framework of the eikonal model. Recall that the unshadowed 
diffractive cross section has a quasibound, 0"°^^ < 2aei ~ 120 mb, valid for very extended 



^We note that ref. obtained 241 mb, i.e., 10% higher than experiment, in their eikonal model. 
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halos, and the actually computed cross section is reduced from this by two factors of two. 
The first reduction, seen in eq. (|T7|), is associated with the fact that the nucleon-target profile 
function blocks a significant fraction of the halo density. The second factor is the shadowing 
of the halo density by the core-target profile function. This shadowing is unavoidable and 
reduces the diffraction cross section by an additional factor of two. 

For the In stripping cross section with an assumed 20% acceptance correction, the eikonal 
model gives — depending on the wave function — between 63% and 83% of the experimental 
result. The cross section increases with larger s-wave contribution in the wave function. But 
we believe that the model that comes closest, s50, is unrealistic on other grounds. Until the 
cross sections are better understood, it will not be possible to use this sensitivity to test the 
s-wave probability in the wave function. 



B. Relative energy spectrum 

In the analysis made in two Breit-Wigner resonances were fitted, at 0.21 ± 0.05 and 
0.62 ± 0.10 MeV, and the relative amounts of s- and p-waves in the wave function were 
extracted. Let us see how this compares with our analysis with the calculated distributions. 
As discussed in Sec. |VT|, we find a peak for the final s-wave at a much lower energy than 
the experimental spectrum shows. However, the finite angular and energy resolution of the 
experimental detectors will inevitably smear out the distribution of the extracted relative 
momentum or energy, shifting the apparent peak to higher energy. This should be taken 
into account in comparing our spectrum with the experiments. 

The authors of ref. give a table of transverse momentum and longitudinal velocity 
resolution widths, and we have folded these widths with the theoretical spectra to compare 
with experiment Q. The folding is done by Monte Carlo sampling. We generate events, 
sampling the theoretical distribution of E in the center-of-mass system of the neutron-core 
system, and assuming that the angular distribution is isotropic there. We add a Gaussian- 
distributed center of mass momentum. As we do not have any cuts in our simulation, this 
center of mass momentum is not important. We transform then from the projectile frame 
to the lab frame, determining the transverse momentum as well as the Lorentz f3 in the lab 
frame. To these we add random errors from a Gaussian distribution with a full width at half 
maximum given by Table 2 in [0. The relative energy (eq. in (3.2) of 0) is reconstructed 
using these values, which are then binned into the same energy intervals as the experiment. 

Fig. O shows the result of this procedure for the s23 model. The experimental resolution 



increases the position of the s-wave peak from 0.08 MeV to 0.15 MeV, which is close to 
the value 0.21 ± 0.05 deduced by the Breit-Wigner fit. However, with our shapes for the 
individual components, the measured peak at 0.21 MeV is a combined effect of both s- and 
p-wave. The individual spectra are strongly overlapping, and require a realistic model of the 
shapes to separate them with confidence. 



^ Ref. [6] also quotes an energy resolution function, (AE ~ E^''^^) but this does not go to a finite 
value near £' = 0, so we prefer to construct the resolution as described. 
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In our treatment, there is no indication of a sharp j9-wave resonance. Our shape has a 
long tail which explains most of the higher energy data. This fact can be demonstrated by 
fitting the experimental data at low energies directly to this shape. One then finds that this 
curve alone is capable of explaining almost all of the higher energy data, leaving no room 
for a resonance there. 

One of our main goals is to find out what we can conclude from this spectrum about the 
s-wave contribution in the ground state and also about the position of the p-wave resonance. 
We investigate this in the following way. The shapes of the distributions are assumed to 
be the same as we found in Sec. |V|, using the parameters of the s23 model. We allow the 
ratio of s to p to vary, as well as the position of the p-wave resonance (that is Eji in eq. 
(^)). Using this as the input to our Monte-Carlo code we calculate the simulated spectrum. 
Minimizing the with the experimental result, we get the most likely result together with 
some range which we still consider likely. Fig. |T3| shows the best fit to the experiment we get 
in this way. It is for an energy of Eji = 0.45 MeV and an s-wave component of about 40%. 
From the variation of "we estimate the uncertainty in the position to be 0.1 MeV and the 
s-wave contribution between 30% and 50%. From Fig. 11 we estimate the s-component in 
the ^^Li wave function to be between 20% and 40% and the most likely value to be about 
30%. 



C. Transverse momentum distribution 

In order to test our final state spectrum, we also make comparisons to a measurement 
of the transverse momentum distribution of neutrons. Such a measurement has recently 
been performed at the same energy p2]. In our analysis we assume that the momentum 



distribution of the neutrons is isotropic and identical to the relative momentum distribution 
of the neutron and the core fragment. The relative energy spectrum can then be transformed 
into a transverse momentum distribution: 

da r°° ,, 1 da 
dki- 



d^k I JO 2'Kmk dE 



In Fig. |1J we compare the result of our "best" model with the experimental result. In 
this figure we made no attempt to include the experimental resolution in the theory curve. 
The agreement between theory and experiment is reasonable especially for medium and 
larger momenta, whereas we have difficulties at the very low momenta. At these momenta 
we expect our approximation of using the relative momenta instead of the neutron momenta 
to be invalid, as center of mass correction should play a role. Also in this analysis we have 
ignored diffraction. 



X. CONCLUSIONS AND OUTLOOK 

In this paper we applied the eikonal theory to nuclear reactions of Borromean nuclei, 
with the object of developing a quantitative tool for interpreting reaction cross sections. 
Because the full theory is quite demanding from a numerical point of view, we examined the 
approximations that are commonly made. In particular, the theory simplifies if one neglects 
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correlations in the ground state wave function, the final state interaction, the distortion 
effects of the profile functions on unstripped particles, or the difference between vectors 
referred to various center-of-mass systems. 

The easiest cross section to interpret is the one-particle stripping, which leaves a particle 
and the core in a final state of low excitation energy. The integrated one-particle stripping 
cross section can be calculated with rather rough approximations. The correlations in the 
ground state play no role except to determine occupation probabilities for the shell orbitals. 
In the differential cross section, the distributions of s- and p- waves are quite distinct, allowing 
the occupation probabilities of the ^^Li halo orbitals to be extracted with 10-20% accuracy. 
Many of the simplifying approximations can be used here without significant error. The 
distortion introduced by the profile functions have practically no effect on the shape of the 
distributions, and a moderate effect on the extracted relative probabilities. However it is 



important to include the final state interaction in the particle-plus-core system ||13[. We 
proposed parameterizations of the s- and p-wave distributions that take both the initial 
halo character and the final state interaction into account. The s-wave distribution, given 
by eq. (pOf), is derived from the asymptotic wave functions. The p-wave distribution, given 



by eqs. (pl| - p2D , is in the form of a Breit-Wigner function with an energy-dependent width. 
The width varies with energy as F ~ E^^'^ at threshold, but the threshold region is very 
narrow due to the extended initial wave function. At higher energies the width grows more 
like r ~ i?^/^, which is the characteristic behavior in the absence of a barrier. 

The other two integrated cross sections are more sensitive to the correlations in the wave 
function. The two-particle stripping cross section is sensitive to the pairing correlation, 
nearly doubling when pairing is included. If the two-nucleon stripping could be measured 
well, it could be used to study this aspect of the wave function. The other integrated cross 
section, the diffractive cross section, is not only sensitive to the eikonal distortions, but we 
found that it changes significantly if the recoil correlation is included in the wave function. 
The differential diffractive cross section is presently beyond our computational powers. 

As a specific example, we have applied the eikonal analysis to the data of ref. 0. The 
total two-neutron removal cross section is larger than our calculation, and also larger than 
the experimental value at a higher energy (800 MeV/n). This is surprising because the 
difference in nucleon-nucleon cross sections at the two energies demands the opposite trend 
in the eikonal model. Our calculated one-neutron stripping and diffractive cross sections are 
only 75% and 50% of the measured values, respectively. The experimental diffractive cross 
section is nearly high enough to be in conflict with the bound given by eq. (|T^). 

In analyzing the differential cross section, we found that both s and p-waves are needed. 
The p-wave cross section does not give a sharp resonant peak, due to the fact that the initial 
wave function is largely outside the centrifugal barrier. The s-wave distribution peaks at a 
very low energy, and the experiment does not show as much of a threshold enhancement as 
we predict. While we are not able to get a perfect flt to the data, our best flt has a p-wave 



resonance near 0.45 MeV, only slightly lower than the measurement of ref. |jT^. The peak in 
the stripping distribution found by ^ is at 0.2 MeV, and is most likely due to the combined 
effect of s- and p-waves. Again, we found that simplifying approximations of no recoil and 
the transparent limit give a rather good account of the one-particle stripping distribution. 
Our analysis with the full eikonal theory gives an s-wave component of between 30% and 
50% in the cross section, which corresponds to an s-wave probability in ^^Li between about 
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20% and 40%. This range also gives a good fit to the transverse momentum distribution of 
neutrons measured by ref. |3^. A number of other experimental |p,17,32|,|6|,|33| and theoretical 
studies |T0| , |3^ have extracted s-wave probabilities of about 50%. An even larger value is 
apparently obtained by ref. [|l^], who report a p-wave probability of 26%. 

In conclusion, the eikonal theory has considerable promise for interpreting the distribu- 
tions and cross sections in nuclear breakup, but it has not yet proved to be a quantitative 
tool for the reaction we studied. 
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TABLES 



TABLE L Parameter of the different models used in comparing the different cross sections. 
The first model ("p89") uses the same potential for s and p wave, whereas the other two use a 
deeper potential for the s-wave (and all other even L waves). The strength of the p-wave potential 
is essentially fixed by the position of the pi/2 resonance. The lowest two entries give the properties 
of the uncorrelated pure s and p-wave models described in the text. 



Model 


Eb 


Vs (MeV) 


ao (fm) 


% Si/2 


% Pl/2 


(fm^) 


(ri - r2)^ (fm^) 




p89 


-0.295 


-35.4 


+L7 


4.5 


89.1 


29.4 


42.8 


18.7 


s23 


-0.295 


-47.5 


-5.6 


23.1 


61.0 


37.7 


45.9 


26.2 


s50 


-0.292 


-51.5 


-90. 


49.9 


33.9 


53.8 


70.1 


36.2 


s 








100 




45.0 


90.0 


22.5 


P 










100 


27.5 


55.0 


13.8 



TABLE II. Integrated cross sections (mb) with different models of the halo nucleus ^^Li. 





uncor. s 


uncor. p 


p89 


s23 


s50 


data 


i 


modified data 


diffraction 


38 


26 


27 


33 


40 


50± 10 


77 


In-stripping 


174 


123 


121 


137 


162 


170 ± 20 


195 


In-str. (transp.) 


182 


129 


134 


155 


182 








2n-stripping 


4 


3 


6 


9 


10 


60 ±20 


8 


2n-removal 


216 


152 


154 


179 


212 


280 ± 30 


280 ± 30 



TABLE III. Integrated cross sections (mb) for He breakup on C. The experimental number 



is from ref. [31|. 





with Sc 


240 MeV/n 
w/o Sc 






with Sc 


800 MeV/ 
w/o Sc 


n 


dif 


32.3 


64.5 


< 2(7ei 


= 119 


42.7 


89.7 


< 2aei = 168 


In-st 


136 


286 






144 


319 




In-st, trans 


170 


409 




= 414 


184 


488 


« 2are = 493 


2n-st 


16.7 


61.6 






19.8 


84.1 




0"-2n 


185 


206 




acxp = 189 ± 14 
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FIGURES 




FIG. 1. Single-particle density of halo neutrons in various models of Li. 
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FIG. 2. Cross sections for nucleon-^^C scattering and for ^Li-^^C scattering. Comparison of 
the results of our profile function with experimental results. 
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FIG. 3. Breakup reaction cross sections for ^^Li-^^C scattering. The left-hand histogram shows 
twice the n-^^C cross sections as required by equahty eq. ([I^) and the bound for the diffractive 



cross section, eq. (18). We compare them with the total cross sections without core shadowing 
("o"'^") in the next histogram. The shadowing reduces the cross section by a factor of ~ 0.4 for 
single-neutron stripping and diffraction, and by a factor of ~ 0.2 for two-nucleon stripping, as show 
by the diagram on the right side. 
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FIG. 4. Impact parameter dependence of the In- and 2n-stripping and the difFractive cross 
sections. Also shown is the square of the profile function of the core-target interaction, S^. 
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FIG. 5. Shadowing factor as a function of target mass number. Shadowing factors are shown 
for 2n-stripping, In-stripping, and diffr active cross sections at an energy of 280 MeV/n. 
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FIG. 6. The ^^Li In-stripping energy distributions in the no-recoil, transparent hmit (solid 
curves) are compared to fits obtained from eq. (|20|) and eqs. (21,2^) (dashed curves). Shown are 
the s- and p-wave components of the s23 wave function. 
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FIG. 7. The ^He In-stripping energy distributions calculated in the no-recoil, transparent 
limit (solid curves) are compared to fits obtained from eq. (2C) and eqs. (21,22|) (dashed curves). 
The s-wave components has been scaled by a factor of 20. 
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FIG. 8. Comparison of the full calculation (including shadowing of the spectator neutron) 
with the transparent limit. Shown are the total In-stripping energy distribution as well as the s- 
and p-wave components of the s23 model. 
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FIG. 9. Comparison of the In-stripping energy distribution of Li with calculations with 
fixed impact parameters of the nucleon-core system (6 and 7fm) and the stripped neutron (1 and 
2 fm). The different spectra have been normalized to the same total value. 
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FIG. 10. Percentage of s- and ^?-wave in the s23 wave function (wfct), compared to the per- 
centages remaining in the neutron-core final state following In-stripping. The different histograms 

show the results of the full calculation and the transparent limit integrated over all impact param- 
eters as well as at specific fixed values of the impact parameter of the nucleon-core system (6 and 
7 fm) and the stripped neutron (1 and 2 fm). 
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FIG. 11. Percentage of s- and p-wave in the various models, comparing the probabiUties in 
the wave functions (wf) to the probabihties in the neutron-core final states {a). 
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FIG. 12. Energy distribution of the s23 wave function. The histogram is the calculated result 
including the experimental resolution. The solid line is the result for all partial waves. Also shown 
are the s- and p-wave components (dashed and dotted lines). The data points are from ref. Q. 
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FIG. 13. "Best" fit of the energy distribution compared to the data of ref. The ratio of s- 
and p-wave and also the position of the p-wave resonance was varied. 
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FIG. 14. Transverse momentum distribution including all partial waves of the best fit model. 
Shown are the results with 30,40, and 50% s-wave component in the final state, corresponding to 
a ^^Li-ground state component between 20 and 40 %. The data points are from ref. 
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